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Abstract. The marked increase in the use of composite and sandwich material systems in 
aerospace, civil, and marine structures leads to the need for integrated Structural Health 
Management systems. A key capability to enable such systems is the real-time reconstruction 
of structural deformations, stresses, and failure criteria that are inferred from in-situ, discrete- 
location strain measurements. This technology is commonly referred to as shape- and stress- 
sensing. Presented herein is a computationally efficient shape- and stress-sensing 
methodology that is ideally suited for applications to laminated composite and sandwich 
structures. The new approach employs the inverse Finite Element Method (iFEM) as a general 
framework and the Refined Zigzag Theory (RZT) as the underlying plate theory. A three- 
node inverse plate finite element is formulated. The element formulation enables robust and 
efficient modeling of plate structures instrumented with strain sensors that have arbitrary 
positions. The methodology leads to a set of linear algebraic equations that are solved 
efficiently for the unknown nodal displacements. These displacements are then used at the 
finite element level to compute full-field strains, stresses, and failure criteria that are in turn 
used to assess structural integrity. Numerical results for multilayered, highly heterogeneous 
laminates demonstrate the unique capability of this new formulation for shape- and stress- 
sensing. 


1 INTRODUCTION 

The inverse problem of shape- and stress-sensing is manifested by reconstruction of 
structural displacements, strains, stresses, and failure criteria using real-time strain 
measurements. Aircraft wings with embedded conformal antennas and those of morphed 
capability require real-time shape sensing to provide feedback for their actuation and control 
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systems. For structural health monitoring, shape- and stress-sensing technologies are the 
enabling capabilities for assessing structural integrity and cost-efficient maintenance. For 
composite and sandwich structures, structural health monitoring using embedded optical-fiber 
networks presents an attractive technology for in-situ strain measurements that give rise to a 
large amount of strain data. Despite their numerous advantages, composite structures may 
experience such modes of failure as delamination and impact damage, and these can affect 
their load carrying capabilities. Thus, the monitoring of structural integrity of multilayered 
composite and sandwich structures is an issue of primary importance. 

Various shape-sensing approaches for plates undergoing bending deformations have been 
explored [1-3]. Bogert et al. [1] examined a modal transformation method, which requires a 
large number of natural vibration modes. Using classical bending assumptions, Jones et al. [2] 
performed the reconstruction of plate deflections by fitting discrete measures of the bending 
curvatures and then integrating the regression curves. Nishio et al. [3] explored shape-sensing 
of thin laminated composite plates by enforcing compatibility between analytic and measured 
strains in a weighted least-squares sense. Both [2] and [3] employed Kirchhoff plate theory, 
thus restricting their methods to homogeneous or nearly homogeneous thin plates. Tessler and 
Spangler [4] proposed a general framework for full-field reconstruction of displacements, 
strains, and stresses, using arbitrary positioned strain sensors on the load-carrying structural 
surfaces. The methodology is based on a least-squares variational principle and accounts for 
the complete set of First-order Shear Deformation Theory (FSDT) modes. They include 
stretching, bending and shear deformations. The variational principle [4] is also well suited 
for finite element approximations. In [5], the authors proposed an inverse Finite Element 
Method (iFEM) based on C°-continuous kinematic approximations. This resulted in an 
efficient three-node inverse shell finite element called iMIN3. 

Although generally regarded as an accurate theory, FSDT may lead to somewhat 
inadequate predictions when applied to relatively thick composite and sandwich structures. 
Fligher-order equivalent-single-layer theories [6] provide improved predictions for such 
structures, specifically for the global response quantities such as deflection and natural 
frequency; nevertheless, even these theories fail to predict through-the-thickness distributions 
of displacements, strains, and stresses with sufficient accuracy. Layer-wise theories [7] 
usually lead to highly accurate response predictions; however, these are obtained at the 
expense of computational efficiency and modeling complexity, especially for multilayered 
structures, since the number of unknowns depends on the number of material layers. The 
recently developed Refined Zigzag Theory (RZT) [8,9] is a good compromise between 
adequate accuracy and computational efficiency. For plate analysis, RZT has seven kinematic 
variables regardless of the number of material layers, just two more than FSDT. This new 
theory is able to model the cross-sectional distortion that is typical of multilayered composite 
and sandwich structures. 

In this paper, the Tessler-Spangler [4,5] iFEM formulation is reformulated to include the 
kinematic assumptions of RZT [8,9]. The new formulation is thus intended for applications 
dealing with multilayered composite and sandwich structures possessing a high degree of 
anisotropy and heterogeneity. The variational principle is then discretized using a C°- 
continuous three-node inverse plate finite element. Numerical results are presented for 
moderately thick sandwich laminates subjected to various boundary and loading conditions. 
Finally, superior stress-sensing capabilities of the present formulation are demonstrated for a 
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select set of challenging material systems. 

2 KINEMATIC ASSUMPTIONS OF THE REFINED ZIGZAG THEORY FOR 
PLATES 

Herein the kinematic assumptions of the Refined Zigzag Theory (RZT) for plates are 
briefly reviewed. In particular, the strain field is formally re-written in order to define the 
strain measures to be used in the iFEM variational formulation (see Sect. 3). 

Consider a plate of thickness 2 h made of N perfectly bonded orthotropic material layers 
(see Figure 1(a)); the superscript ( k ) denotes the kth layer. The plate is referred to a Cartesian 
coordinate system (x 15 x 2 , z ) where (x, , x 2 ) are the in-plane coordinates and z is the 

thickness coordinate that ranges from -h to h , with z = 0 identifying the mid-plane and z {J) 
identifying the yth interface (see Figure 1(b)). 

The displacement field of RZT for plates is [9] 

u[ k \x l ,x 2 ,z) = u(x l ,x 2 ) + z 0 l (x l ,x 2 ) + <fi < l k) (z) y/ l (x l ,x 2 ) (1) 

u ( 2 k) (x l ,x 2 ,z ) = v(Xj,x 2 ) + z 6 1 (x ] ,x 1 ) + (j) < 2 ) (z) ^ 2 (Xj,X 2 ) 

U_ (Xj ,X 2 ,z) = w(Xj , x 2 ) 

where u[ k) and u 2 k) are the in-plane displacements and u z is the transverse displacement. 
RZT has seven kinematic variables, u = [«, v, w, O v 0 2 , i//, , y / 2 ] 7 . u , v, and w are the uniform 
displacement components along the x, , x 2 , and z -axis respectively; 6 { and 0 2 are the 
average rotations of the transverse normal around the positive x 2 -axis and the negative x, - 
axis, respectively; and y/ a (a- 1,2) are the amplitudes of the zigzag contributions to the in- 
plane displacement in the x a -directions (see Figure 1(a)). The zigzag terms <f> (k) i// a (« = 1,2) 

in Eq. (1) describe the C° -continuous cross-sectional distortions that are typical of multilayer 
laminates. The zigzag functions, <t> (k) {z) , have units of length and are piecewise linear, C°- 
continuous functions of the thickness coordinate and of the transverse shear moduli of the 
laminate layers. Refer to [9] for the detailed derivation of the zigzag functions. 
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z (0) 
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x l (x 2 ) 


Figure 1 : Plate notation (a) and layer notation for a three-layer laminate (b) 
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The in-plane components of the strain field are given as 
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represent the membrane, bending and zigzag strain measures, respectively. Note that, whereas 
the membrane and bending strain measures are constant with respect to the thickness 
coordinate, the zigzag strain measures have the “zigzag”, z-dependent distributions. 

The transverse shear strains are given by 

fi ( * )= *7i + ( 1 + A w yi (4) 

A? = Vi + { 1 + 02 k) )v' 2 


where J3 ( a k) = ( a = 1,2) and 

SU U l/j -Mi l W ,2+^2-^2j 


with rj a ( a = 1,2) denoting the transverse-shear strain measures of RZT. 


(5) 


3 INVERSE FINITE ELEMENT METHOD BASED ON RZT 

In this section, the iFEM variational formulation based on the RZT kinematics is 
presented. The formulation enables reconstruction of the deformed shape of composite and 
sandwich structures from in situ strain measurements. 

The general framework is that of finite element approximations. Thus, a discretization of 
the structure with plate elements is introduced, in which the element kinematic variables are 
interpolated by a set of suitable shape functions, 

u(x!,x 2 ) = u e =N(x 1 ,x,)q e (6) 

where N(x,,x 2 ) denotes the shape functions and q the nodal degrees-of-freedom. 
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3.1 Error functional 


Following the iFEM methodology [4,5], the displacement solution is obtained through the 
minimization of an error functional, which is defined as the least-square error between the 
analytic strain measures (Eqs. (3), (5)) and their measured values, known at discrete locations 
from in-situ strain measurements. For a single element, the error functional is given as 


<r( u e ) = |e(u e )-e £ 2 + k(u e ) 
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where the superscript “e” is used to denote the measured values; for the membrane, bending 
and zigzag contributions, the squared norms that appear in Eq. (7) have the form 
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where n is the number of locations P = (x l ,x 2 ) i , where the strain measures, ef , k' ; and nT (/) 
are evaluated from strain-sensor measurements. The zigzag strain measures, m[ (/) , are 
evaluated at the /th interface using embedded strain-sensors. The transverse shear strain 
measures, g(u), cannot be obtained experimentally. Thus, the transverse shear term in Eq. 

(7) is given by the L 2 norm 



( 9 ) 


where A e denotes the element area. In Eq. (8), a is a positive valued, small (compared to 
unity) weighting coefficient; refer to the numerical studies in Sect. 4 in which the value of 
a = 10 -5 is used. 

Considering Eqs. (3) and (5), while invoking the kinematic field interpolations given by 
Eq. (6), the strain measures can be expressed in terms of the nodal degrees-of- freedom, q e , as 

e ( ue ) =B e(*i> x 2)q e > K ( ue ) = B K( x i’ x 2)q" (io) 

m (<) (u e ,z) = B m (x 1 ,x 2 ,z)q‘ > , g (u e ) = B g (x, , x 2 ) q e 

where the matrices B e , B K , B m and B„ contain the derivatives of the shape junctions 

N(x I 5 x 2 ). 

Substitution of Eq. (10) into Eqs.(8) and (9), and minimization of the functional cp e with 
respect to the unknown nodal degrees-of-freedom, q e , yields the element matrix equation 
a q = b , where the matrix a e is given as 
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a? = - Z [ B c ( X 1 ’ X 2 ), B e (*1 » *2 ), ] + Z [ B K ( X 1 » ' *2 ), B K (*1 > *2 ), ] 
^ 1=1 ^ 
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whereas the right-hand-side vector is 
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Taking into account appropriate coordinate transformations, the element contributions are 
assembled into a global system of equations. Upon enforcement of problem-dependent 
displacement boundary conditions that prevent rigid-body motion, the inverse-problem 
equations take on the form 

Aq = b (13) 


where A is a well-conditioned square matrix. 

Equation (13) is well suited for real time applications. This is because A is inverted only 
once (assuming small displacements, the strain-sensor locations, [(x p x 2 ),,±h\ and 
[(Xj ,x 2 ) j ,z U) ] , remain unchanged.) On the other hand, the vector b needs to be updated at 
each strain-data acquisition increment. Thus, the displacement solution is efficiently 
computed by the vector-matrix multiplication, q = A b . The strains and stresses are readily 
computed for each element using strain-displacement and constitutive relations of RZT. 
Furthermore, strains and stresses can be used to construct failure criteria. 


3.2 Evaluation of the experimental strain measures from strain-sensor data 

The present method is especially aimed at structures with embedded strain sensors. In 
particular, a strain-sensor configuration that measures the in-plane strains is considered on the 
top and bottom surfaces and at one interface (refer to Eq. (2)). Considering that the zigzag 
contributions to the in-plane strains vanish on the top and bottom surfaces [9], the expressions 
for the experimental membrane and bending strain measures are derived by evaluating Eqs. 
(2) at the discrete locations P t = (x 15 x 2 ) ; . , z = ±h. 
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where jX , e 22 , ;/, + 2 } and \s u ,s 11 ,y n \ denote the in-plane strains measured respectively at 
[(Xj ,x 2 ) ( .,+/z] and [(x p x 2 ) i ,-h\. Measuring the same strain components at the /th interface, 
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{ £ n’ £ 22 ’Tn} > an( l again using Eq. (2), the zigzag strain measures are determined as 
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Using this strain-sensor configuration, the zigzag strain measures, mf ( , are evaluated at the 
/th interface only. 

4 THREE-NODE PLATE INVERSE ELEMENT BASED ON RZT 

In this section, a brief description of a three-node inverse plate element is presented. The 
element has seven degrees-of-freedom at each node, u k ,v k ,w k ,6 xk ,6 lk ,y/ lk ,y/ 2k , where 
k = 1, 2, 3 is an index ranging over the three nodes. Considering that only first derivatives of 
the kinematic variables appear in the functional, Eq. (7), the shape functions are required to 
satisfy C°-continuity. For this purpose, the so-called anisoparametric interpolations are used 
(e.g., refer to Tessler and Hughes [12] and Versino et al. [11]). 

The in-plane displacements, bending rotations and zigzag amplitudes are interpolated using 
linear shape functions. The functions are defined in terms of the area-parametric coordinates 
L k (k = 1,..,3), 

3 

x{x l ,x 2 )=Y J L k{ x u x 2 )Zk^ x={u,v,e l ,e 2 ,i// l ,i// 2 ) ( 16 ) 

k = 1 

whereas a quadratic interpolation is used for the deflection 

3 3 

w(x 1? x 2 ) = ^L k w k + y\L ]k (6 lk — yS\k) + L 2k ( d 2k y^ 2 k ) J (17) 

k = 1 k = 1 

where the quadratic shape functions, L Xk and L lk , are given by the expressions 

L U = y(M,, - a m L , ) L lk = ^ (b m L, - b,L m ) (18) 

a k={ X Un~ X u) b k =( X 2l- X 2m) 

the subscripts being given by the cyclic permutation of k = 1-3 , / = 2, 3, 1 and m = 3, 1,2. 
The anisoparametric interpolations, Eqs. (16) and (17), ensure that (a) truly thin plates can be 
modeled without any stiffening due to shear locking, and (b) the resulting element, herein 
referred to as iRZT3, has the same number of degrees-of-freedom as a standard linear- 
interpolation element. 

5 NUMERICAL EXAMPLES 

A symmetric three-layer sandwich plate, with carbon-epoxy face-sheets and a PVC core is 
analyzed, under the action of different loadings and boundary conditions (see Table 1 for the 
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mechanical material properties). The thickness of each face-sheet is 10% of the total laminate 
thickness. In the following example problems, square plates are considered having edge 
length a and span-to-thickness ratios, p = a/2h = \0. Regular mesh patterns are used having the 
same number of elements along the plate edges, n e . Figure 2 shows the plate geometry, the 
coordinate system, and a 2 x 2 discretization (n e = 2). 

Two example problems are analyzed: (1) a simply supported plate subjected to a transverse 
bi-sinusoidal pressure p(x l ,x 2 ) = p 0 sin(7rxja)sin(7rx 2 /a) , and (2) a cantilever plate 
subjected to a uniform transverse pressure p(x l ,x 2 ) = p 0 . In lieu of the actual experimental 
strain measurements, an exact elasticity solution by Pagano [13] is used for the simply 
supported plate. For the cantilever plate, an analytic RZT-based direct solution is used, [9]. 
In both cases, the strain sensor configuration consists of three strain-sensor rosettes located at 
the element centroid on top and bottom surfaces and at Z( i> (see Figures 1 and 2). 

Table 1: Mechanical properties of orthotropic (C) and isotropic (P) materials. The Young’s moduli and the shear 

moduli are expressed in GPa. 
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Figure 2: Inverse plate finite element discretization, 2x2 mesh (n e = 2), and strain-sensor locations. 

In Figures 3-5, the in-plane displacement, in-plane normal stress, and transverse shear 
stress are normalized as 

“i = u l x(l0 4 D n /p o a 4 ), <t„ = < T n x(\/p 0 p 2 ), f 2z = t 2z x (\j p 0 p) (20) 

where Z), , is the bending stiffness coefficient. 

Figures 3-4 depict results for the simply supported plate where the displacements and 
stresses obtained using the present element formulation (labeled as iRZT3) are compared to 
Pagano’s exact elasticity solution. In Figure 3(b), for comparison purposes results are also 
depicted for the predecessor inverse element “iMIN3” [5], based on FSDT. When using n e = 
12 (Figure 3(a)), the present iRZT3 model accurately predicts the maximum deflection (error 
within 0.27%) whereas iMIN3 leads to an underestimation by 70% (not shown). Figures 3(b), 
4(a) and 4(b) show that through-the-thickness distributions of in-plane displacement, in-plane 


8 


Priscilla Cerracchio, Marco Gherlone, Marco Di Sciuva and Alexander Tessler. 


normal stress and transverse shear stress are accurately recovered using iRZT3 whereas 
iMIN3 in unable to model the zigzag shape of the in-plane displacement (Figure 3(b)) typical 
of a sandwich-like stacking sequence. 



(a) (b) 

Figure 3: Simply supported plate, a/2h=\0: (a) Percent error of maximum deflection, e w = 100(w/ w Pasa "° -1) , 
vs. the number of elements along the plate edge, n \ and (b) Through-the-thickness distribution of in-plane 
displacement, u t (0,a/2) , for the discretization n = 12. 




Figure 4: Simply supported plate, a/2h=\0: through-the-thickness distributions of (a) In-plane stress, 
a n (a/2,a/2 ) , and (b) Transverse shear stress, T 2 .(0,a/2) , for the discretization n= 12. 
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In Figure 5, results are depicted for the cantilever plate, where the iFEM predictions are 
compared with the corresponding RZT analytic solutions of the direct problem. Figure 5(a) 
depicts rapid convergence of the maximum deflection error, whereas Figure 5(b) shows a 
through-the-thickness distribution of the normal stress, cf, , (a/4, a/2) . 
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(a) (b) 

Figure 5: Cantilever plate, a/2h=\0: (a) Percent error of maximum deflection, e w = 100 (w/ w' ar —1) , vs. the 
number of elements along the plate edge, n\ and (b) Through-the-thickness distribution of in-plane stress, 

cf,, (a/4, a/2) , for the discretization n= 20. 


6 CONCLUSIONS 

Real-time reconstruction of structural displacements and stresses from in-situ discrete- 
location strain measurements, herein referred to as shape- and stress-sensing, is an inverse 
problem that has important implications for monitoring of structural integrity, as well as for 
the actuation and control of smart structures. A new variational formulation for shape- and 
stress-sensing of laminated composite and sandwich plates has been presented. The approach 
makes use of an inverse Finite Element Method (iFEM) that was previously developed on the 
basis of First-order Shear Deformation Theory (FSDT). The iFEM is based on the 
minimization of a least-square error functional in which analytic and measured strains are 
enforced in the least-square sense. Within the present formulation, the kinematic assumptions 
of the Refined Zigzag Theory (RZT) are used and require that the measured strains be 
available along three surfaces through the laminate thickness (instead of only two surfaces for 
the FSDT-based formulation.) The RZT has previously been shown to be especially well 
suited for the modeling of laminated composite and sandwich plates. The error functional is 
discretized using C°-continuous interpolations of the displacement field, yielding an efficient 
three-node inverse plate finite element that has seven kinematic degrees-of-freedom at each 
node. From the reconstructed element displacements, strains and stresses at every material 
point of the structure are then computed using RZT’s strain-displacement and constitutive 
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relations. Numerical results for moderately thick laminated sandwich plates undergoing 
elasto-static deformations demonstrated superior full-field predictions for the displacements, 
strains, and stresses. As expected, the present iFEM-RZT formulation results in more accurate 
predictions than those based on iFEM-FSDT, and this is especially evident in highly 
heterogeneous and sandwich laminates. 
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